LISA: Locally Interacting Sequential Adsorption 

_j Anton Muratov* Sergei Zuyev* 

O March 1, 2013 

(N 



(N 






X 



Abstract 



We study a class of dynamically constructed point processes in which at every step a 
new point (particle) is added to the current configuration with a distribution depend- 
ing on the local structure around a uniformly chosen particle. This class covers, in 
particular, generalised Polya urn scheme, Dubbins-Freedman random measures and 
r^ cooperative sequential adsorption models studied previously. Specifically, we address 

^ models where the distribution of a newly added particle is determined by the distance 

C^ to the closest particle from the chosen one. We address boundedness of the processes 

I I and convergence properties of the corresponding sample measure. We show that in 

general the limiting measure is random when exists and that this is the case for a wide 
'"^ class of almost surely bounded processes. 
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1 Introduction 

A model of sequentially constructed point process that inspired this paper was 
presented to one of the authors (SZ) by Richard W. R. Darling as a way to 
describe a certain population dynamics. His original model is described as 
follows. Start with a fixed finite configuration X = {xi, . . . ,a;„j,} of no points 
in a plane. Call them particles. Choose one of these particles uniformly at 
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random. This particle, say ^, is thought of as a 'parent' of a new particle 
that will be added to the current configuration according to the following rule. 
Consider k closest to ^ particles a;i(^), . . . ,Xk{^), where fc > 3 is a parameter 
of the model, and fit a 2-variate Normal distribution centred in ^ to these. Let 
V be the corresponding estimate of the covariance matrix. Then sample a new 
particle from this estimated law: Xng+i ~ A/'(^, V). Once this is done, we have 
a configuration of rig + 1 particles and we repeat the procedure again: choose 
randomly a particle among all uq + I particles now present, estimate the Normal 
law from the closest to it k particles and add a new particle sampled from this 
law, etc. 

A realisation of the model based on 20 initial particles after 10 thousand 
steps is shown in the upper-left plot in Figure [Tl 

One can note the following characteristic features of the construction. Since 
the parent particle is chosen uniformly, there is a greater chance that this parent 
will be chosen in the area densely populated by the particles. Moreover, in 
these dense areas the distance between the particles tends to be small, so the 
newly added particle also tends to lie close to the parent point. So as the 
construction progresses, it tends to reinforce the dense areas of particles which 
kind of 'adsorb' new particles. This is clearly seen in Figure [T] where the 
configuration is shown after 10, 20, 30 and 40 thousand steps. Each of newly 
added 10 thousand particles are shown in dark emphasising their trend to follow 
higher density areas of the previously existing (grey) points. Note that although 
the configuration of existing particles plays a crucial role in the construction, 
only k closest particles to the chosen one actually contribute to the distribution 
of the added point. In this sense the interaction is local, hence the name we have 
chosen for this process: Locally Interacting Sequential Adsorption or LISA, for 
short. 

Another feature concerns the geometry of the cloud of particles. When the 
parent particle is chosen inside a circular cloud, its closest neighbours tend 
to be homogeneously spread around it. This produces more or less isotropic 
Normal density for a new point wich adds to a round cloud making it even more 
isotropic. In contrast, when a boundary particle is chosen as a parent or when 
it lies in a stretched cloud, the density will also be skewed in the corresponding 
direction. So in the long run round clouds tend to stay round, but time to time 
'shootouts' from their boundary happen which then tend to produce filamentary 
arrangements. It also happens due to randomness that even if a parent in such a 
filament is chosen, it can still produce a particle well outside the main direction, 
and this would then become a centre of another circular cloud. 

There is a range of questions arising immediately: will the particles be al- 
ways confined to a bounded region or will the diameter of the cloud will increase 
indefinitely? Will eventually particles be present in any compact set of a pos- 
itive area or will there be gaps never filled by the process? If we supply all n 
the particles present at the current step with masses 1/n we get a probability 
sample measure f„. Is there a limit in appropriate sense of the sequence of 
these measures? Is this limit measure when exists random or is it non-random? 



And what about finer properties of this hmiting measure, hke the Hausdorff 
dimension of its support? 

Surely, the two-variate Normal distribution is just one of possible choices of 
the distribution governing addition of a new point. And all the above questions 
can be asked for any other distribution: for instance, to provoke shootouts one 
would try some heavy-tailed distribution for the distance from its centre. We, 
however, want to keep the main essence of the local interaction of the model 
above requiring that the new particle distribution scales appropriately when the 
configuration becomes denser. This will bring us to the notion of a stopping set 
described in details in the next section. 

The structure of the paper is the following. In the next Section [2] we fix 
the notation used throughout and give formal description of the class of locally 
interacting sequential processes we are dealing with, the Darling's model being 
one particular case of these. Other cases include such seemingly different models 
as Dirichlet measures, Dubbins-Freedman's random distribution functions and 
cooperative sequential adsorption. Section [3] demonstrates on a simple exam- 
ple that the limiting distribution of particles, if exists, is generally a random 
measure, this particular example leads to Dubbins-Freedman construction of a 
random distribution function. Section [4] addresses boundedness issue and show 
that under rather mild conditions the cloud of points has almost surely finite 
diameter. Finally, Section [5] studies convergence of sample measures and shows 
that in models with an a. s. finite diameter such a limiting measure exists in a 
weak sense almost surely. LISA processes constitute a very large class of models 
with different properties, so we conclude by outlining extensions, relations to 
other models and open problems which are abound. 

2 Preliminaries and Model Description 

In order to define a locally interacting sequential adsorption process, we need 
a few components. First, the phase space W, where the particles live, and an 
initial configuration Ar„^ == {xi, . . . , x„(,} of particles in it which is a parameter 
of the model. Although a generalisation is immediate, we assume in this paper 
that M^ is a subset of Euclidean space K'^. It is often convenient to treat a 
particle as a unit mass measure so that a collection of particles is a counting 
measure on the Borel subsets of W. 

As already alluded in Introduction, the local interaction, thought of as a 
dependence of the distribution of the newly added particle on the local config- 
uration of particles around its parent, can be described in terms of a stopping 
set which is the next component to be defined now. 

Let M denote a set of Radon measures on the Borel sets B of M'* and Af C Ai 
be the set of counting cr-finite measures on B. For a closed set G & B, let J"g 
be the cr-algebra of subsets of M generated by the sets {/x € A4 : fi{B n G) < 
t}, B e B, t > and let J" = VbJ'g, where G runs through any countable 
system of bounded closed sets generating B. The system {J^g} is a filtration^ 
because it possesses the following properties: 



1. Monotonicity: Tc C Tq whenever G' C G; 

2. Continuity from above: Tq = ^n^^Gn fo^' ^-^y sequence of closed nested 
sets: Gi 3 6*2 3 . . . such that n„G„ — G. 

A random measure (resp., a point process) is a measurable mapping from some 
probability space to [A^, J^] (resp., to [Af,J-]). A realisation of a point process 
is called a configuration (of particles). 

Denote by F the ensemble of all closed sets of M.'^ and by S the smallest 
(T-algebra containing the sets {G g F : G n X 7^ 0} for all compact sets K. A 
random closed set is a measurable mapping from a probability space to [F,S]. 
We will be working with the canonical space for the point processes when dealing 
with random sets so they become a measurable functions of point configurations. 

A stopping set is a random closed set S : [Af, {J^g}, P] '^ Wi 2] such that the 
event {S C G} is J^g measurable for any G € F. The corresponding stopping a- 
algehra Fs consists of events E ^ F such that E\^{S C G} G J-q for any G S F. 
A stopping set is a generalisation of the classical notion of a stopping (or Markov) 
time: likewise a random process' trajectories stopped at the Markov time, the 
geometry of a stopping set is determined by the configuration of particle inside 
it and on its boundary and does not depend on the particles outside of the 
stopping set. 

For more details on stopping sets in M"^, see [9] and Appendix in [Tj covering 
also more general phase spaces. 

Returning to the construction of LISA, secondly, for any point x & W and 
all finite configurations X with n > no points there is defined a stopping set 
Sx{X \ {x}) (by definition, X \ {x} = X ii x ^ X). In other words, if X' is 
another configuration such that X' n Sx{X \ {x}) — X {^ Sx{X \ {x}), then 
necessarily Sx{X' \ {x}) = Sx{X \ {x}). From now on, to ease the notation, 
we will simply write S{x^X) of just S^ when no confusion occurs instead of 
Sx[X\{x}). 

Finally, for every stopping set Sx with the corresponding stopping cr-algebra 
J-S^ there is defined a random variable Cs^ , such that its distribution is de- 
fined only by the geometry of the stopping set Sx and the particles it contains. 
In other words, Sx can be viewed as a parameter of this distribution, or if 
there are other natural parameters of this distribution, they are necessarily 
J\g^ -measurable. Typically, for our purposes Sx and Cs^ are defined to be shift 
invariant and scale homogeneous, so that 

S{x,X)^x^ S{{),X -x) Cs(x,x)^x + Cs(Q,x-x) (1) 

5(0, aX) ^ a 5(0, X) Cs(o,aX) = a Cs(o,x) (2) 

for any positive a, configuration X and x (z W {= denotes equality in distribu- 
tion). In R. Darling's model described in the previous section, the stopping set 
S{x, X) is the smallest closed ball centred at a; e A" containing k nearest neigh- 
bour particles of A \ {x} to x. The covariance matrix V estimated from these 
particles (with or without x itself) defines (s^ as having Multivariate Normal 



distribution MVN(x, V) centred at x. Since only the particles contained in Sx 
are used to estimate V^ V is J-5^ -measurable. 

Having these necessary components at hand, we define a dynamical proce- 
dure by which new particles are sequentially added to the existing configuration 
one by one. Let {x„} be a sequence of independent random variables, where Xn 
is uniformly distributed on the discrete set {1,2, ... ,n}. Given current config- 
uration Xn = {xi, . . . ,Xn} oi n > uq particles, a new particle Xn+i distributed 
S'S Cs{x ,x„) is added to the configuration. In other words, a particle of X„ is 
uniformly chosen (so it is a particle with index Xn), and then a new particle is 
added according to the distribution defined by its stopping set. 

We now give examples of models which are constructed this way. 

Example 1. Let W = [0,1], no ~ 1 and Xn^ — {0}. The stopping set Sx 
is the segment from x to the next particle to the right (or to 1 if there is no 
such particle). Formally, Sx — S{x,X) = [x,x + d'^{x,X)], where d'^{x,X) = 
min{a;' — X : x' G X\{2:}U{1}, x' > x}. Finally, Cs^ is a uniformly distributed 
point on Sx. 

In the example above all the added particles belong to W = [0, 1] by con- 
struction. In the next model the particles' range grows indefinitely, but as we 
show in the next section, all the particles will be confined to an almost sure 
bounded (but random) set. 

Example 2. Let VK = M and Xng be some set of no > 2 particles. The stopping 
set Sx = [x — d{x,X),x + d{x,X)], where d{x,X) = m.in{\x — x'\ : x' G X\{x}} 
is the distance to the closest to x particle of X. As in the previous model, Cs^ 
is uniformly distributed in Sx- 

Exam,ple 3. This is one-dimensional variant of the model described in Introduc- 
tion. Here W, Xng , Sx and d{x, X) are as in the previous example. But ^s^ 
is Normally distributed with mean x and standard deviation ad{x, X) for some 
a>0. 

Exam,ple 4. More generally, let W = W^, S{x, X) he the closed ball centred in x 
with radius d{x, X) and (s^ — x + d{x, X)tp with a given random vector ijj e M'' 
whose distribution does not depend on anything. In the previous two examples, 
d = 1 and ip is uniformly distributed in [—1, 1] in Example^ or -0 has normal 
A/'(0,a^) distribution in Example [s] 

In the next two examples, the distribution of the nth new particle to be 
added does not depend on the index variable Xn , but rather on the whole current 
configuration of the particles. 

Exam,ple 5. Let W he some measurable space and /i be some given probability 
measure on its measurable subsets. Define S{x, X) to be the whole W for 
all X and X. Random variable Cs^ equals x with probability 1 — 1/n and 
otherwise a random variable with distribution /i with probability 1/n, where 
n is the cardinality of X. Surely, the parameter n of the distribution of Qs^ 
is J"s^ = Fw measurable. Such defined LISA process describes the Blackwell- 
MacQueen construction which generalises the Polya urn scheme. It weakly 
converges to a Dirichlet random measure in the limit, see [3]. 



Example 6. Let W be some compact subset of K" and {Pn} be a given sequence 
of positive numbers. Fix also a positive parameter R called the interaction ra- 
dius. Define Sx to be a closed ball B{x,R) centred at x with radius R and 
Cs(x,x) to be the random variable with the density proportional to the func- 

tion f{x) — f{x,X) — Y[k=o f^ri{x,x), where n{x,X) is the number of particles 
from X belonging to B{x,R). The corresponding LISA process then defines 
the so-called cooperative sequential adsorption (CSA) model, see [51 15] and the 
references therein. 

When the stopping set Sx is allowed to be the whole W, we are basically 
in the situation when the distribution of the added particle depends on the 
whole current configuration. Such construction may include just about any 
dynamically constructed processes and is too general to be treated in a unified 
manner. So to stay in the "locally interacting" framework, we will concentrate 
in this paper only on LISA processes where the distribution of ^5^ depends only 
on the distance (i(a;, X) from x to the (properly defined) closest particle among 
X \ {x}, i. e. on Examples Il]-l4] The original Darling's model which inspired this 
investigation does not fall into this framework (unless fc = 1 in ID case) and 
its detailed analysis is still a hard open problem. But even the models we do 
analyse here exhibit fascinating and different behaviours. These concern, first 
of all, randomness of the limiting distribution, boundedness of its support and 
its dimension. 



3 Random limiting distribution of LISA 

This section demonstrates that, in general, the limiting distribution of LISA pro- 
cesses is non-degenerate. We show on Example [T] that the sample distributions 
functions F„(i) = n~^ 122=1 ^Xk<t, t S [0, 1] of the first n particles converge to 
a random distribution function on [0, 1] arising in the Dubbins-Freedman con- 
struction, see [4]. This fact has already been noted in [71 Sec. 5.2], but included 
here for a didactic purpose. 

Recall the Dubbins-Freedman construction of a random measure with sup- 
port on [0, 1]. A realisation of the cumulative distribution function of such a 
measure is produced by the following sequential procedure. Let ui — {(pi^ij^i) 
be two independent uniformly distributed in [0, 1] random variables, or, equiva- 
lently, Ui ~ Unif([0, 1]^). The vertical and the horizontal lines passing through 
this point divide the square [0,1]^ into four rectangles. The distribution func- 
tion being constructed is deemed to pass through the points (0, 0), {(pi, ipi) and 
(1, 1). Since the c.d.f. is a non-decreasing function, it must be contained in the 
rectangles: Ii^i = [0,<pi] x [0,-!/'i] and /i^2 = ["i^i,!] x [-01,1] lying along the 
'main' diagonal from bottom left to top right. Namely, the c.d.f. passes through 
a uniformly generated point uii ^ Unif(/i 1) in the first rectangle and through 
a uniformly generated point ui_2 ^ Unif(/i_2) in the second. Again, both these 
points divide the corresponding rectangles /i^i and /i_2 into 4 rectangles each, 
the diagonal ones containing the c.d.f. In each 4 of these diagonal rectangles of 
level 3 random uniform points are selected which the c.d.f. is deemed to pass 



through, etc. Thus the c.d.f. is defined on a everywhere dense set in [0, 1] and the 
values in ah other points are defined as the hmits. Thus one obtains a continuous 
increasing curve which is a random element on the space of independent uniform 
random variables indexed by a binary tree: Ui,ui^i, ui 2, ui 1 1, ui 1 2, Wi,2,i, ■ • ■ • 

Consider now the first particle X2 generated in the construction described 
in Example [1] It has Unif[0, 1] distribution. Now an analogy with Polya urns 
can be drawn in the following manner. Paint the particle cci = black and 
the particle X2 white. The next generated particle X3 will be black or white 
according to whether it is the black xi or it is the white X2 selected on stage 
3, i. e. X2 = 1 or X2 — 2. Then the procedure repeats with the colours of the 
particles being inherited from their 'parent' particles. So the number of black 
and white particles has the same distribution as the number of black and white 
balls in a Polya urn with starting configuration of one black and one white 
balls. But all black particles are lying to the left of X2 and all white are to 
the right. So the proportion of the black particles is the proportion of particles 
with coordinates less than X2 which equals the proportion 2/2 of black balls in 
the urn scheme which is Beta(l, 1), or equivalently, the uniform distribution on 
[0, 1]. Thus the limiting c.d.f. passes through the points {x2, 2/2) having the same 
distribution as ui in the Dubbins-Freedman construction. 

Conditioning now on the value of the second 'daughter' particle x' of Xi = 0, 
the proportion of all the particles to the left of it conforms to Unif [0, 2/2] • Indeed, 
just ignore all the white particles in the construction and distinguish among all 
'black' particles the ones which are really black in [0, x') and 'dark grey' which 
lie in [x',X2)- So the c.d.f. passes through the point distributed as ui_i above. 
Iterating to other segments, we conclude the demonstration of the equivalence. 

4 Boundedness of LISA processes 

Next natural question to be addressed is whether the limiting distribution of 
particles in LISA, when it exists, has an a. s. bounded support. This is trivially 
true for Example [l} but it is not that evident in other examples. Notice, that 
since the series 1/n diverges, each particle will be chosen infinitely many times 
as a parent point. Thus the rightmost particle present at stage n, for instance, 
in Example [2] will eventually be chosen and with probability 1/2 will produce 
a particle yet more to the right. So the support is growing, but will it stay 
compact nevertheless? 

4.1 Boundedness in Example [2| 

Recall Example [2] Let W — R, X„g = X2 = {0, 1}. It is convenient to slightly 
reformulate the rule by which new particles are added: 



f-n+l 



^n+ldnVn+l , (3) 



where {£«} is a sequence of i. i. d. random variables equal to ±1 with probability 
1/2, {r]n} are i. i. d. uniformly distributed in [0,1] and given x — x^^, dn = 



d{x,Xn) G J^s{x.x„) for the stopping set S{x,Xn) — [x — d{x,Xn),x + d{x,Xn)]- 

Theorem 1. Denote 

rUn — niinja; : x £ Xn} , 
Mn — maxja; : x e Xn} 

for the LISA model in Example [2| Then 

P{— oo < liniinf m„ < limsupAf„ < +00} = 1 . (4) 

Proof. Wc only prove that hnisupM„ is finite a. s. A proof of finitencss of 
liminf rn„ is similar. Introduce i>m — time of the rn-th jump of the process 
{Mn} as follows: 

i/Q == no , 

Vjn+i = minjfc : Mk > M,,^}, m = 0, 1, 2, . . . 

Now consider imbedded process M^ — M^^^ . From (Is]) we obtain 

M„i = ^xum.-! + ^''n.-i£i'„,^i',„> '71 = 1, 2, . . . (5) 

The distribution of rjn is concentrated on (0,1], thus the maximum can only 
have m-th jump at time n if Xn — ^m- That means, on the (i^nj— l)-th step 
the {m — l)-th maximum is chosen, implying that x^^ -^ — Mm-i- Moreover, 
e^^-i must be equal to 1 in order for a positive jump to happen. Thus ([5]) is 
reduced to 

Mm = Mm-1 + dy^rju^, TO = 1, 2, 3, . . . (6) 

Notice that d^^ — d^^{M,n-i) is less or equal than M^-i ~ A/„i_2, and for 
that expression we have 

Mm-i - Mm-2 = d^^-iVu„,-i, m = 2,3, ... 
By induction, one can get 

rn 

Mm-Mm-l < n^'^'^' 
fc=l 

so, since Mq = 1, we come to a bound 

m I 

M™<1 + ^[]77,,,TO = 1,2,3,... (7) 

1 = 1 k=l 

Notice now that {f]i,^}k>i and {rik}k>i are equally distributed, since by the 
definition of Vm, n^ ^ vi < V2 < vz < . ■ ■ , and {77^} are independent of {vm}- 
The sequence {Mm} is monotonely increasing, therefore it has an a. s. limit 
Moo, although possibly infinite. However, since E?7„ < 1, 

E lim Mm — lim EAf^ < 00, 

TTl— >oo ?n— foo 

in particular. Moo < cjo a. s. thus finishing the proof. D 



Remark 1. Same proof with minor tweaks works for any initial configuration 
Xno = {a;i,a;2, . . . ,Xno} and even for rj in the scheme ([8]) distributed with an 
arbitrary law F concentrated on (0,1] with P{r]n < 1) > (to guarantee a 
convergence of the sum below) . Estimate (l7| turns into the following: 

rrt I 

M,n < dno (Mo) (Mo + ^Y[vu,),m^ 1,2,3,... 

1=1 fe=i 

Here Mq and d„g(Afo) are constants dependent on Xng only. 
Remark 2. Dubbins-Freedman random distribution functions F(t) arising in 
Example [I] are almost surely continuous, but also each point t G [0, 1] is almost 
surely a point of growth of the distribution function F{t), i. e. for any e > there 
exist t' G {t-e,t) and t" G (t,i + e) such that F{t') < F{t) < F{t"). In contrast, 
Example [2] provides random distribution functions which are continuous, but 
also contain constant regions, i. e. the corresponding limiting measure does not 
have connected support. To see this, observe that with positive probability 
there happen to be a configuration of points generated by the algorithm of 
Example [2] where there are 2 pairs of points: xi < X2 < yi < 2/2 each pair 
consisting of closely situated points separated by a relatively large void, i. e. 
X2 - xi,y2 - yi are small but (j/i - X2)/(x2 - Xi) and {yi - X2)/{y2 - 2/i) are 
large. The construction of new points scales with distance, so that the evolution 
of the initially present pair of points at distance S has the same distribution as 
the evolution of two initial points at the distance 1 scaled by factor 5. Therefore 
according to just proven Proposition [l] with a positive probability the maximum 
of all the offsprings of the pair a;i,a::2 (affected only by the distance X2 — X2) 
will be strictly smaller than the minimum of all the offspring of the pair yi , 2/2 
(based only on 2/2 — 2/1) so that there will be a void somewhere between X2 and 
2/1 not filled with any points. 

4.2 Boundedness in Example [4[ 

Important feature of the model considered in the previous section is that the 
particles cannot jump over each other and thus the influence of a new added 
particles can be effectively controlled. This is now longer the case in Example |4] 
where the farthest particle can be potentially generated by any parent point. In 
this subsection we present sufficient conditions, under which the more general 
d-dimensional LISA process from Example [4] is bounded a. s. 

Let W = R"^, d G N. Initial configuration is given by Xq — {xi, a;2, . . . , Xng}. 
New particles are added according to the rule: 



I'n+l 



^Xr, + c^«(a^x,.) V'n+i, n = no, no + 1, . . . (8) 



As before, x„ are independent, distributed uniformly on {0, 1, . . . ,n}, (i„(x) is 
the distance from x to Xn\{x}, ipn are i. i. d. random variables with a given dis- 
tribution which may now have a non-compact support. Next, we set r]„ = \\'4'n\\ 
and denote C = E771. Also, put ijn = r/nAl = min{77,i, 1} with corresponding 
C = E,r}i. 



Lemma 1. Let {??«}, {'fin} be i. i. d. sequences of non-negative random variables, 
independent between themselves. Put fin — rjn A 1, 0n = ^n nr=i ^i- ^^^ ^^^'^ 
C = E771 and E(^i be Rnite. Put C = E771. Assume C < 1. Define 

n 

Then Yco = max Yn < 00 a.s. Moreover, 

n 

Eyoo<--^(1 + E^i) (9) 

Proof. Note that Yn arc monotone and thus converge to some (possibly infinite) 
hmit. However, one can show that Yn are bounded in C^: 

71 n n i—1 

Ey„ - E Y 0,(^, + 1) < ^E0,(<^, + 1) = ^Er?i [] E%(E(^, + 1) = 

2—1 i—1 i—1 :? — 1 

Hence yoo has finite expectation with a correspondent bound. D 

Theorem 2. If C + C < 1 then sup |^„| < 00 a.s. Moreover, 

n 

Esup|g„|<ylo+ ^ V ^ , 

n 1 — G — G 

for some constants Aq, Dq depending only on the initial conhguration Xn„. 

Proof. First of all, organize {a;„} in a tree in the following natural way: for every 
point Xig from X„„ denote all the points Xn such that Xn~i — io as {xi„i-^}°^^i, 
in the order of appearance. We will further say that Xi^i-^ are the children of 
Xig. Then all the children of Xi^i-^ we denote by {x 1^1,^12} i°^i, and so on. 

Let i>{io . . .ik) denote the (random) time of appearance of Xi^.i^.. Let also 

tpzo-.tk — '0,.(io. ..«*.)' Vto-tk — Vi^(to...ik)^ c^^o...^fc(•) — d^(io...tk)-ii-) ■ Observe 
that {ipio...ik}kAo,...,ikeN and hence {rii„,„ik}kAo.--;-i.kefi are i. i. d. families. 
Fix io < no for now. We have in our new notation: 

Estimate (iiQij^(xij,), ii = 1,2,... Recall that it is a distance to the closest 
neighbour, hence it can not be larger than distance to the points that already 



10 



exist for sure at the moment of Xj^jj 's appearance, that is, its mother Xi„ and 
all of the older sisters Xi„i,Xig2, ■ ■ ■ , 2:jo,«i-i- We can write: 

digiyxig) < ang[Xig) =: dig 

dio2{x^„) < d,„i{xig) A \x,„i - x,J < di,,i{x,„) A |di„i(a;,JV'»oi| 

< di„i{xig){l A \'^p^ol\) < di„fi,„i 

dio3{Xi„) < di„2iXio) A \Xi„2 - X.i„\ < dio2{Xi„) A \dig2{Xio)i>ig2\ 

< d^g2iXig)il A |V'»o2|) < digfl,gifl,g2 



n-1 
"io*i(-''*o) — dia YY ^*oJ 



(11) 



Introduce 0io...ij^ := 77io...i^,_ji^ H Vio-.-ik^ij- Using ( 1 1 ) we can estimate the 



right part of (10): 



n-i 



V \dioil{Xio)'4'ioil\ < dio y IV-ioill n '"'*0J "" ^'0 V'''"'! n ^*0J 

Apply Lemmalllwith 0„ = ^ionj <^n = 0, n = 1, 2, . . . to see that L^° is a proper 
random variable with E L^° < — ^ . 
Now estimate the second generation. 

Y \Xi0ili2 ~ Xiol S Y I I'^^oil ~ Xiol ~r \f \Xioi1i2 ~ ^Jj^j^ | I 
il,i2 il \ 12 / 

— V I I'^'oii ~ Xiol ~r Y "ioiii2 (.-^ion J'7ioni2 I l-'-^J 



Note that 



22 — 1 



dioiii2\Xioii) — "io«l 1 i-^ioil ) XI ^*o*iJ — 

J=l 



■^ioii "^io 



12-1 

in 

i=i 



Vioilj 
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and therefore we continue ( 12 ) 




< V 1^*0*1 - 2;»ol 1 + \J nio'ii'i2 n ''^o^iJ 



j=i 



<rf»oV(^^on (l + V^*o*i*^)) 






Note that {V^Joiii2}i"=i i^ ^^ i.i.d. sequence, distributed Hke V^ and inde- 

pendent of {^ioJi}ir=i- Using Lemma 1 again, we obtain V2 < 00 a. s. , and 
moreover, 

- " 1-c Vi-c. 

Repeating that argument for n = 3, 4, . . . we obtain a monotone sequence 
{Ln}?^=i: where digU° is an a. s. bound for elements from n-th generation 
of descendants of the point Xi„, n = 1,2, .. . Since {LJi"} is monotone and is 
bounded in >C^: 

f-[\i-c) l-C-C 

there exists an a. s. limit L^ := lim U° with EL^ < — 77— a- We finish the 
proof by recalhng arbitrariness of iq: 

■nil I I -n I \ A ^0-^0^ 

Esup|x„|< sup Ixjj+Esup sup Ix^o...^^ - x^J < Aq H -^ — — 

n>l iQ<nQ ri iQ,...,?^! 1 — Cy — G 

Here Aq = sup l^ij, i^o = sup dn^ixi^) D 

As an illustration, the ID Darling's model in Example [3] is bounded if a < 
0.8239 which is a value obtained numerically. 

Remark 3. The bound in ^ and therefore the condition C + C < 1 may be 
improved if one could find "nicer" conditions sufficient for 

EiJ = Emax(?7fcf)i?72 • . .?yfc-i) < 1 

for i. i. d. non-negative rjk- We demonstrate that for a particular distribution of 
rjk- If we introduce C^ = log?7fc, with the distribution function F(t), then the 
calculation shows that G{t) — P{\ogH < t) must satisfy the integral equation: 

Cit\ = / ^-oo G(^ - ^) ^^W + (^W " FmG{t), t > 
^^ \/!^G(i-z)dn^), t<0. ^ ^ 
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Let 77fc have the foUowing distribution: 

'O, 2/<0, 

P(%<y) = <((i-p)/, o<2/<i, 

for some a>l,/3>0, 0<p<l. In that case, Cfc = — log r/h has c.d.f. 

^(,) f(l-p)e^ ^<0, 



and one can directly obtain the solution for (13 1, 



leading to the following distribution for H: 

(0, x<0, 

i(^2:-" + l)-i-^, x>l. 

It is not possible to find an analytic form for E iJ in that generality, but if we fix 
/3 to be, say, 1/4, then we can find the regions where EH < 1 and E 77 + E 77 < 1 
numerically: 




60 SO 



So, as we see, the condition E?7 + Ef/ < 1 of Theorem |2] is far from being tight 

for boundedness. 



5 Properties of the limiting measure 

In previous sections we addressed the boundedness of the series of point con- 
figurations X„. In this section we study the limiting sample measure and its 
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properties. We are still working with most general model of Example [4] in the 
phase space W = K'', d > 1 and initial configuration X„(, — {xi, . . . , Xng}- New 
particles are added according to the rule: 

x„+i = x^^ +d„ix^^)i;n, n = n+l,n+2,. . . 

so that Xn+i — Xn U {xn}. Again, Xn are independent, uniformly distributed 
over {l,2,...,n}, ipn are i. i. d. random variables distributed according to a 
given probability measure ^, dn{x) is, as before, a distance from x to its closest 
neighbour in X„ \ {a;}. Denote by 



M« = - V Mn,fc = - V ^J.{d';,'^{ ■ - Xk)) 

k=l k=l 

the distribution of a;„+i, and denote by 

1 " 

the empirical measure of the process X„ after n steps. 

We will start with a short lemma providing some insight on the behaviour 
of dn{x). 

Lemma 2. Assume that \jXn = X^o is a.s. bounded and P(|V'i| < 1) > 0. 

n 

Then 

lim max dn{xk) = 

n— >oo l</c<n 

Proof. First, notice that for every A:, lim dn{xk) = 0. This follows from 

n— >oo 

P(|V'i| < 1) > and X]P(Xn = fc) = oo, i.e. every point is going to be 
picked up an infinite number of times and infinite number of times its dn{xk) is 
going to shrink. 

Next, assume the contrary. Let 

lim sup max dn{xk) > e > 

ji l<k<n 

Pick {kj}, {rij} such that 

dn, {xk, ) > £, J e N 

Since dn{xk) monotonely tends to zero as n goes to infinity, we can assume all kj 
to be different and moreover, dk{xk) > dn{xk)- That means, in particular, 
that 

\xki -Xk^\ > e, j eN, i < j, 

i. e. Xoo can't be covered with a finite number of balls of radius s — a contra- 
diction with the a. s. boundedness. D 
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Lemma 3. Assume that [J X„ = Xqo is a. s. bounded. Assume that one of the 

n 

weak Umits lim /^„, lim z^„ exists a.s. and is equal to /i*. Then the other one 

n— ^oo n— ^oo 

exists a.s. and is equal to fi* , too. 

Proof. Let p be the Levy-Prokhorov distance between probability distributions. 
We will use the following property: 

py^Q-kFk.^akGkj < niax(^l,^afej niaxp(Ffc,Gfc) (14) 

Taking it into account, one can write: 

1 " 1 " \ 

'^^ 6{- - Xfc), - V" ^„,fe < max p {5{- - Xk),p.n.k) < 

71 ^ — ' / l<k<n 

fc=l / ^ ^ 

p((5(-),/i( max d^^(a;fc)-)) ^ 0, n ^ cxi, 

l<k<n 




smce 



lim max dn{xk) — 0, 

n— ^oo \<k<n 



and for every p, - probability measure on M, p{an-) — )■ 5, whenever a„ — > cx). D 

As we have noted, Example [T] is equivalent to the Dubbins-Freeman con- 
struction, so the limiting measure p* exists. The Hausdorff dimension of its 
support, in a slightly more general setting, was found in [S], which is equal to 
1/2 here. We are going to use the technique from [7] to show that a limiting 
measure exists in Example [2] however, it is still an open question, how to cal- 
culate the Hausdorff dimension of its support and whether a limit exists at all 
in Example [3j 

Now, we have W^ = M, -0„ - Unif(-l,l). Let {xn,k}l=i,n = 1,2,... 
denote the rearrangement of the elements of X„ in ascending order. Then the 
complement of X„ consists of n-l-1 intervals, In,j'- 

In,o = (-oo,a;n,i) 

J-n,k ^ [Xn,kj Xn^k+l)j fc ^ i, Z, . . . , n i. 

Lemma 4 (cf.Lemma 2.1 in [T). Let < b < a be integers, F C (0,1,2,..., a), 
cardF = b. Then 

Mn(|J Ia,j) -^ Z{a.F) 
j€F 

Here Z(^a.F) is independent of {xj}i<j<a, and its law has density 

-t^-\l-tf- 



r(a) .b-1/1 ,\a-fc-l 



r(&)r(a-&) 

forb = b- i(I(0 eF) + l{a e F)). 
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Proof. Follows from the result on the generalised Polya's urn scheme: probabil- 
ity for a new point to appear in J„j- is 2;^ if it is /„.o or /„^„ and - otherwise. D 

Corollary 1. 

lim sup lim ^m{In 1) = 

Theorem 3. Almost surely exists ^* - a probability measure such that 

Proof. By Lemmaplwe can define almost surely an increasing function G{x) on 
Xoo: G{x) :— lim /i„i((— 00, x)). Note that since almost surely 

—00 < iniXoo < swpXrx, < +00, 

thus sup(a; : G{x) = 0) > —00 and inf(x : G{x) = 1) < +00. Moreover, by 
corollary, 

sup(G'(a;) : x < y) = inf(G(x) : x > y), 

therefore, almost surely G{x) is a continuous cumulative distribution function 
for some probability measure /i*, which is easily shown to be a weak limit of 

{Mn}. □ 

We will now prove a couple of facts about the third model. First, let us make 
the following observation. Denote by d* the maximal spacing of a configuration 

d* = max dn{xi) 

l<i<n 



Theorem 4. If Vn — XI ^xi Is an empirical measure of the process Xn on the 

i=l 

n-th step, then the Levy-Prokhorov distance between the empirical measures 
for the two consecutive configuration is given by the following expression 

p{vn, v„+i) = min(n"\ max(n"^, d*J) . 

We will now present some bound for the decrease of the maximal spacing. 
The technique we use is essentially due to [7j. 

Introduce ip{t) = 1 - {Erf + E77*) and a = sup {^}. 

t>0,V5(t)<l 

Theorem 5. If C + C < 1, then limsupn'^d* < 00 a. s. 

Proof First step of the proof is to estimate dn{xk) from above with a certain 
well-behaving construction. Introduce a triangular array {A„jj} I < k < n, n > 
1 as follows. 

A„o,fc = dng{xk),k = 1,. . .,no, 
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For n > no 'we put 

^n.kVn, iffc = n + l, 
An.k otherwise. 

In other words, at each step pick one "diameter" and replace it by two diameters, 
scaled with the realisations of ry and 77. That corresponds to the new point being 
added to the configuration at the n + 1-th step with its initial distance, and it's 
mother point's distance being scaled correspondently. 

It is not very hard to see that for any n > nQ,k < n, one has a bound 

d„{xk) < A„^fc. (15) 

Now, prove that the configuration {A„ j,} behaves nicely. Pick a positive t so 
that Et]* +E7/* < 1. Such t exists, because of the starting conditions. Consider 
the quantity X:fe<„A^,fe- 

Y. ^n+l,k = Y. ^i» - ^'Vx„ + ^n,X„ (^« +Vn)^J2 ^nM ~ ^n,xM^^) 
k'Cn+l k<n k<n 

Then if J^n is the sigma-algebra generated by the evolution of X„ up to time n, 
one has 



k<.n-\-l k<.n 



and so Y{i<j<n{^- ^) ^Yjk<n^n,k^ together with sigma-algebra J"„, is a 
positive martingale, which has an a. s. finite limit. However, 



Mt) 



l<3<n •' l<j<n ■' -^ 

has a positive limit as well, because < (p{t) < 1. Therefore 

fc<n 



converges to a finite limit as n -H- cx). That, together with (15), finishes the 



proof. n 

We conclude our paper with noting that all of our models can imbedded 
into the continuous time in a natural way: each particle produces children 
independently according to a Poisson process with intensity 1, and new points' 
distribution depends only on the local configuration around the mother point at 
the time of birth. The embedded processes of particle births is equivalent to the 
original LISA process. Indeed, when n particles are present at any given time t, 
because of the lack of memory of the exponential distribution, the next one to 
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give birth is uniformly distributed among them. The model in Example [T] then 
becomes the so-called fragmentation or stick breaking process, see, e. g.,0 and 
the references therein. In continuous time version of LISA, the total number 
of particles is a continuous time Galton- Watson process of pure birth. In that 
case, one can also obtain an asymptotic bound for the maximal spacing of the 
point configuration. 

Theorem 6. If Xt is the continuous time version of the third model, d^ is the 
maximal spacing of Xt , then 

limsupe~(i-'^-^)*Edt < oo. 

£— foo 

The proof is essentially representing {^t}tg[o,T] ^-s a trajectory- wise limit 
as n goes to infinity of a discrete-time process {X(m);}i<;<^„7-j in which each 
particle at each step gives birth to a new one with probability Pm = —, and 
then observing that if we implement {^i,/c}i<;^i„y| similarly to the previous 

proof, so that A/^^ is a bound for di{xk), then one can write out the following 
distributional inequality 

Ar+i ^ (1 - Ii)A: + MA^fii A A^ryO < (1 - Ii)A: + h{K*if^i + Kl\) 

where A^^ is an independent copy of A^. Therefore, one can argue that 

EA^,„^j<(l-(l-C-CW)L™^J 

which in the limit brings us to 

EA:<e-(i-'5-C')T^ 

6 Conclusion and open problems 

We have defined a wide class of dynamically constructed point processes where 
new particles are added randomly with their distribution depending on a local 
neighbourhood of a randomly uniformly selected particle. Exact notion of local- 
ity is based on stopping sets methodology. In this paper we considered only the 
case of models where this set is the ball centred in the selected particle with the 
radius equal the distance to the closest particle. Obvious generalisation is to 
consider the ball to the k-ih. nearest neighbour, as it is the case in the original 
Darling's model described in Introduction. But even for the case fc = 2 we were 
not able to find ways to control the spread of the particles' cloud. So all the 
questions of boundedness and/or existence of a limiting measures remain largely 
open. 

Boundedness results we established in Sectionl4]are only sufficient conditions. 
In fact, we were not able to prove unbounded behaviour in any model. Our 
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hypothesis is that if the distribution of the variable ||V'|| in Example [4] is heavy- 
tailed (the shootout distance distribution before scaling), this should produce 
clouds of particles which spread indefinitely. 

As already mentioned in the previous Section, existence of a limiting measure 
is also an open question for all the models but two Considers there. We tend to 
think that boundedness of a model should suffice for a weak limit to exist. 
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Figure 1: Sequence of 10, 20, 30 and 40 thousand generated particles. Newly 
added ten thousand particles are in dark, previously existing - in grey (partially 
covered), initial particles are contoured void circles. 
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